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ABSTRACT 


In order to design a controller for a linear physical plant, the transfer function of the plant 
must be obtained. A system consisting of Texas Instruments TMS320C30® digital signal 
processor, an analog interface board and a personal computer is used as a test bench platform 
for real-time system identification. The recursive least squares algorithm is implemented on 
the TMS320C30® digital signal processor to determine the coefficients of the plant transfer 
function in real-time using input-output data from the plant. The program for identification 
can be called from the Matlab” environment. Three system identification examples are 
presented. 


I. INTRODUCTION 


A system which facilitates rapid identification of the transfer function of a control system 
plant is a useful tool in a control laboratory. A system which consists of a Texas Instruments 
TMS320C30° digital signal processor (DSP), a custom-built analog interface board (AIB) 
and a personal computer has been developed for rapid identification of the parameters of a 
linear plant transfer function. The input and output data from the plant are collected and 
processed in real-time using a recursive least squares algorithm to determine the coefficients 
of the plant transfer function. The Bode plots of the identified transfer function and the 
actual transfer function are compared. 


An interface program on the personal computer makes it convenient to run the identification 
program from the Matlab” environment. The system is described in section II, the 
identification algorithm and software are presented in section III. Results on the 
identification of the transfer functions of a passive first order low pass filter, and a passive 
second order low pass filter, and a second order dc servo system are presented in section IV, 
conclusions are given in section V. 
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Il. SYSTEM DESCRIPTION 


The system for the identification of the parameters of a plant transfer function consists of Texas 
Instruments TMS320C30° DSP evaluation board, a custom-built analog interface board with two 
10-bit A/D converters and two 12-bit D/A converters. The TMS320C30° DSP evaluation board 
equipped with 16K of memory occupies a slot in the personal computer. The TMS320C30® can 
perform 33 million multiply-accumulate operations per second. The personal computer is used 
for developing code, downloading code to the TMS320C30° processor, for collecting the results 
and displaying the transfer function parameters. The identification program is run from the 
Matlab® environment. 


A block diagram of the system set-up for plant transfer function identification is shown in Figure 
1. 


DSP 
EVALUATION 


BOARD 


Figure 1 Block Diagram for System Identification 


The plant input data u(t) and output data y(t) are digitized by the AIB board and sent to the 
TMS320C30® for processing. 


Il. IDENTIFICATION ALGORITHM AND SOFTWARE 


Suppose the continuous-time plant to be identified is represented by a transfer function, 


Y(s) | N(s) 
U(s) D(s) 


G,(s)= (1) 


where N(s) and D(s) are polynomials. Selecting a sampling interval T and assuming that a zero- 
order hold precedes the plant, an equivalent discrete time transfer function of the plant is, 


Y(z)_ bz 't+K+b,z" 
U(z) 1-a,z'-K-a,z” 


G,(z)= (2). 
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The objective is the determination of a, and b, from the plant input-output sequence u(k) and y(k) 
using recursive least squares algorithm implemented on the TMS320C30° processor [1]. The 
least squares algorithm is given in the following steps [4]. Let the estimate of the coefficients of 
the plant transfer function be denoted by a vector, 


1. Select a, y, and N where a=y=weighting factors=1 for least squares, and N samples. 

2. Select initial values for P(N) and i N) were P(N) is an error estimation and i N) are the 
present estimated system coefficients. 

3. Collect y(0),K ,y(N) ,the outputs, u(0),K ,w(N), the inputs, and form ®7(N +1), the arma 
model. 

4. LetkeNn. 


1 Pls DR we ey) 
a new gain correction factor 


5. Let rer Mates af + 
Y a id 


matrix. 
6. Collect a new y(k+1) output and u(k+1) input. 
7. Gk+1) — Gk) + Lk+ (yk +1) -®7(k+ )@(k)) the next estimate of system 


coefficients. 
1 
8. P(k+l= ma —L(k+I@"(k+ 1)] P(x) a new error estimate. 


9. Formanew ®(k +2) using y(k +1) and u(k +1). 
10.Letk=k+1. 
11. Go to step 5. 
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The recursive least squares algorithm in step 1-11 was programmed in C programming language 
on the PC. A C compiler was used to translate the C program into assembly language, and the 
machine code of the assembled program was downloaded to the TMS320C30® processor for 
execution [2]. A program was written to run the identification program from the Matlab” 
environment [5-6]. After the computation of the transfer function parameters, the coefficients are 
displayed on the screen of the PC. A flow chart showing the execution of the recursive least 
squares algorithm is shown in Figure 2. 


Enter System Parameters, i.e. Order 


Recursions, and Initial Conditions 


Download Parameters to DSP 


Compute Plant Coefficients for 


Number of Recursions 


Upload Plant Coefficients 
from DSP to PC 


Display Coefficients 


Figure 2 Flow Chart Showing Execution of Recursive Least Squares Algorithm 


The computer programs for the identification may be obtained from the authors. 
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IV. APPLICATION EXAMPLES 


Three transfer function identification examples are presented in this section. The first example is 


a first order passive low pass filter shown in Figure 3 


R= 100kQ 


u(t) = 16Vp-p 10Hz Square Wave 


Figure 3 Passive Low Pass Filter 
The transfer function is, 


45.45 


H = —__ 
= Sas 


(3). 


Using a sampling interval T = 2 x 10° sec, and a zero-order hold gives the discrete equivalent of 
(3) as, 


009z" 


(2) =F 975 


(4). 


Using a 10Hz square wave of 16Vp-p amplitude as input signal u(t), the identification algorithm 
produced an estimated transfer function of, 


008427! 
A(z) = 
M2=T ogi © 


after 200 recursions. Bode plots which compare the ideal transfer function and the estimated 
transfer function as the number of recursions is varied are shown in Figure 4. 
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Ideal 

10 Recursion 
~ 20 Paecursins 
~ 200 Reoursans 


Magnitude (dB) 


10° 19° 
Frequency (radianssec) 


Phase (deg) 


10° 10° 
Frequency (racians!sec} 


Figure 4 Comparison of Bode Plots of H,(z)and #,(z) 


The next example is a second order low pass shown in Figure 5. 


u(t) = 8V-10Hz Sauare Wave Biko R2=30.06kO 


C1=0.225uF C2=0.231uf 


Figure 5 Second Order Low Pass Filter 


The transfer function is, 


21471110" 


6 
s’ + 4.40957 x 1075+ 2.14711 x10 (6) 


H,(s)= 
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Using a sampling interval T = 1 X 10” sec, and a zero-order hold gives the discrete equivalent of 
(6) as, 


H(z) 0.009301 1z' + 0.0083055z~ 
Zz — 
1-1.62609z' + 0.643422? 


C1): 


Using a 10Hz square wave of 8Vp-p amplitude as input signal u(t), the identification algorithm 
produced an estimated transfer function, 


0.05216z ' + 0.0089013z 


AB(z)- 
2(2) = 71966142 +033040722 


(8) 


after 500 recursions. Bode plots which compare the ideal transfer function and the estimated 
transfer function as the number of recursions is varied are shown in Figure 6. 


and Order System fs=500 Samples par gecond! Pls (4aradsac 6 P3=44radipor 
TT TErer T TT TTT T TT TTT T 


' keual 
‘ —— 10 Recursons 
|} —— 800 Recursion 


Magnitude |) 


Dreecpaa recy [ndet ce | 


102 
Frequency [reais] 


Figure 6 Comparison of Bode Plots of H,(z) and #,(z) 
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The final plant whose transfer function is to be identified is a 2" order servo system whose 
transfer function 1s, 


= 133152 
S° + 26.6885 + 1328.9656 


H,(s) = (9). 


Using a sampling interval T = 1 x 10° sec, and a zero-order hold gives the discrete equivalent of 

(9) as, 

— 0.6598 x 10° z' + -0.6540 x 10° 2” 
1-197242' +0.97372z” 


Hz (10). 


Using a 1Hz square wave of 1 V amplitude as input signal u(t), the identification algorithm 
produced an estimated transfer function, 


0.00677z' — 0.032165z~ 
1—0.610122z7' — 03612482” 


A(z) = (11) 


after 500 recursions. Bode plots which compare the ideal transfer function and the estimated 
transfer function of the number of recursions is varied are shown in Figure 7. 
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Figure 7 Comparison of Bode Plots of H,(z) and F(z) 
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The identification of the 1“ order transfer function H(z) in (4) produced F(z) with parameters 
very close to those of H,(z). Figure 4 shows that the magnitude and phase plots of H,(z) and 
F(z) are very close after 200 recursions. The identification of the 2" order low pass transfer 
function H,(z) in (7) produced an F(z) with coefficients which were not very close to those in 
(7), however, the magnitude and phase plots in Figure 6 show that B,(z) is a fairly good 
approximation of H,(z). The identification of the transfer function of the 2™ order servo system 
was not very successful. The coefficients of the identified transfer function F(z) were very far 
off from those of the ideal transfer function H,(z) in (10). The recursive least squares algorithm 


was unable to identify the resonant frequency of the servo system. Further investigation is 
required to determine the cause of the discrepancy between the ideal and estimated transfer 
functions. 


V. CONCLUSIONS 


A digital signal processor based system for the real time identification of the plant transfer 
function using the Texas Instruments TMS320C30° processor has been developed. The 
TMS320C30° processor computes the coefficients of the plant transfer function from input- 
output data using a recursive least squares algorithm. The computed transfer function 
coefficients are displayed on the screen of the PC. The system is useful for rapid real-time 
identification of system transfer functions in a control laboratory. Further refinements of the 
algorithm are needed to make the system a more accurate identification tool. 
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